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Abstract. The earth's imperfectly 
elastic resoonse to body and loading tidal 
forces is discussed using complex Love 
numbers and complex mass loading coeffi- 
cients. Exact analytical expressions have 
been derived relating the energy dissipa- 
tion within an inhomogeneous, comnressible 
solid earth to the surface values of these 
complex characteristic numbers, thus rela- 
ting the global dissipation function Q to 
the phase shifts in the potential, gravity, 
tilt, strain and displacement tides. 

Integration of a global ocean tidal mo- 
del shows that energy dissipated in the so- 
lid earth due to ocean loading is at least 
10 % of that dissipated in the body tide; 
however both body and loading tide together 
do not account for more than a few percent 
of the astronomically observed dissination. 

The commonly used relation Q =— tan^P 

where ip is the observed phase lag only app- 
lies v/hen selfgravitation and hydrostatic 
prestress are ignored. It, therefore, is not 
applicable to the earth and in fact there is 
no unique relation between the global Q and 
the tidal phase shifts, this relation being 
very dependent on the distribution of Q with 
depth. Determinations of the global Q from 
satellite observations may be in error by 
70 %, and calculations on the basis of seis- 
mic Q-models predict phase shifts in the 
gravity tide of only a few thousands of a 
degree in place of the currently predicted 
tenths . 

Unlike in the body tide case, dissipation 
in the loading tide is sensitive to proper- 
ties of the asthenosphere, and phase shifts 
in the M_ loading tides in displacement, gra- 
vity and'^tilt may be as high as several de- 
grees for loads near ocean ridges and sub- 
duction zones. 

Rate of Tidal Dissipation in the 
Solid Earth 

Dissipation of body and loading ti- 
dal energy within the solid earth may 
be determined from its complex Love 
numbers and complex mass loading co- 
efficients, respectively [Zschau, 1979 a, 
b] . Such the expressions for the dissi- 
pated energy turn out to be fairly simple, 
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for the loading tide, where Ae is the 
energy dissipated during one cycle of har- 
monic loading, K andH* , are the imagi- 
nary parts of the surfSce £ove numbers and 
mass load coefficients, respectively, ijj is 
the amplitude of the body force potentiSl, 
ijj* is the amplitude of the surface load 
potential, R is the earth's radius and G is 
the gravitational constant. n describes 
the degree of the expansion into spherical 
harmonics. The integration is taken over 
the surface of the earth. The expressions 
above have been obtained without approxima- 
ting the real earth by an incompressible 
and homogeneous one as was necessary in for- 
mer calculation, for instance by Munk & 

Mac Donald [i960] . 

Using these formulas, and assuming the 
mantle Q structure LM.S as obtained from 
the observation of the earth's free oscilla- 
tion [see Smith, 1972] , a body tide solid 
dissipation rate of , ,, 

3.19 X lO"^ erg/s , 

i.e. about 1 % of the astronomically obser- 
ved dissipation rate, has been obtained. The 
corresponding computation for the loading 
tide dissipation rate within the earth's 
crust and mantle, carried out on the basis 
of a global M 2 ocean tide model [Hender- 
shott, 1972] , yields a minimum value of 

roughly 10 % of the body tide solid dissipa- 
tion rate, i.e. ,, 

3 . 21 X lO"^ erg/s . 

This value has been determined from the low 
degree harmonics of the ocean tide distribu- 
tion up to n=25, and, therefore, does not 
represent the high amplitude tides in the 
shelf areas. The latter may contribute 
significantly to the total dissipation, 
because the dissipated energy is propor- 
tional to the square of the marine tidal 
amplitude. On the other hand, the elastic 
strain energy stored in the mantle as well 
as the energy dissipated in the mantle de- 
creases with increasing degree n of the 
spherical harmonic loading for n> 5 as 
may be seen from Fig. 1. This suggests 
that solid earth dissipation in the shelf 
areas does not change the total dissiua- 
tion rate drastically. Anyway, both body 
tide - and loading tide dissipation to- 
gether most probably do not account for 
more than a few percent of the astronomi- 
cally observed dissipation rate. 
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Fig. 1. Loadina tide shear strain energy 
stored and dissipated within the solid 
earth [from Zschau, 1979b] . n is the 
degree of the spherical harmonic expan- 
sion of iiendershott ' s global ocean 
tide model [iiendershott, 1972 ] . Results 
are valid for the Gutenberg Earth v;ith 
the free oscillation Q structure LMS for 
the mantle and the surface v;ave 0 struc- 
ture for the crust. 


Body Tide Phase Shifts and the Earth's 
Dissipation Function Q 


In Zschau [1979b] the earth's global 
dissipation function Q has been related to 
the phase shifts in the potential, gravi- 
ty, tilt, strain and displacement body 
tides. The following expressions valid for 
constant Q distribution within the earth 
have been obtained: 


Tangential tan^) = 

displacement: 

Potential va- tan^p = 
riation due to 
tidal defor- 
mation : ' 

Radial dis- tan*P = 
placement : 

Surface areal tanip = 
strain : 


L/1 = - 0.897 Q 
KA = - 0.612 Q 


H/h = - 0.555 Q 


H-3L 

h-31 


- 0.318 Q 


Tilt 


gravity va- 
riation 


tan^= = + 0*217 o“^ 

tan^P= = " 0.0508 Q 


(h,H) , (k,K) , and (1,L) complex Love 

numbers; the index n has been omitted here. 



CONVENTIONAL 

CORRECT 

OBSERVED 

l/tonfli. 1 l 

Qgje 


DELAY ANGLES 

H— 2 K 

(CONSTANT Q) (LOW 0 ZONE) 

RECORDING 0.1° 

79 

2B 

5.5 

GRAVIMETERS 0.2° 

39 

14 

<3 

SATELLITE 

l/tona|, 



OBSERVATIONS 0.5° 

60 

36 

34 


TABLE 1. Conventional and correct global 
Q values (0 .) corresponding to tidal 
delay anglei'^of gravity and satellite 
observations . 


The above equations shov? that the commonly 
used relation ^ 

o" = - tan^P (3) 

where ip is the observed phase shift is not 
applicable in the earth tide case. For the 
tidal gravity variation this had already 
been pointed out by Slichter [1960] who 
instead of (3) introduced the formula 


q"^ = - ^ tantp (4) 

where 6 is the gravimetric 5 “factor, 
and is the gravity phase shift. 

With 5 “ 1.16 this yields 

tanip = - 0.138 q"^ (5) 


which is still too large by more than a 
factor of 2 as is obvious from the compa- 
rison of this formula v;ith the correspon- 
ding one given above. One may show that 
Slichter 's formula is equivalent to 
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which does not correspond to the basic de- 
finition of as the strain energy dissi- 
pated during one cycle of loading over 2Tt 
the peak energy stored in the system. 

The general expression for the earth's glo- 
bal dissipation function Q in terms of its 
complex Love numbers is, however. 



Kh-Hx (1+k) 
h* (1+k-h) 


[Zsch au. 


1979b] (7) 


This equation is valid for an incompressib- 
le body with homogeneous density. None of 
the above phase shifts is equivalent to this 
expression, hence in general 

q“^ * - tanlp (8) 

As shown in Zschau [l979b] the reason for 
this is that in the case of tidal deforma- 
tions, sel fgravitation and hydrostatic pre- 
stress cannot be ignored. 

Besides the fact that (3) is not appli- 
cable in the tidal case, there, further- 
more, is no unique relation between the glo- 
bal Q and the tidal phase shifts, this re- 
lation being very dependent on the distri- 
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bution of Q with depth. For instance, a Born [1976J deduced the phase angle of 

low Q zone in the upper mantle such as given the potential bulge due to the body 
in model LMS [see Smith, 1972] may alter the tide of Mars. They related this nhase 
ratio between observed tidal phase shifts and angle to the global Q of Mars by the 


the global (average) Q up to a factor of 5. 

Table 1 gives some examples for the 
errors involved when' not taking account 
of the above aspects: Let the phase delay 
of the body tide gravity be 0.1° with re- 
spect to the external forces. Proceedlncr 
in the conventional manner, i.e. using 
Slichter's formula with 5 = 1.16 gives 
the wrong global body tide Q = 79. In 
the case of constant Q within the 
mantle, we find the correct value to 
be Q = 28. If we assume a low Q asthe- 
nosphere, i.e. let the real Q distribu- 
tion within the earth differ by only a 
constant factor from the free oscilla- 
tion Q model LMS, we find that the global 
Q has to be chosen as low as 5.5 to 
correspond to the gravity phase delay of 
0.1°. This is less than 7 % of the value 
Q = 79, obtained by the conventional me- 
thod. For the Seune reason it turns out 
that the body tide gravity phase delay 
due to friction within the solid earth 
amounts more likely to a few thousandth 
of a degree than to a few tenth of a 
degree as expected so far. We, therefore, 
suggest that the average 01 gravity phase 
shift of - 0.2° as observed for Europe 
may not be attributed to imperfect 
elasticity in the mantle like it is pro- 
posed by Melchior et al. Il976], but 
rather to the indirect effect of the 0. 
tide in the oceans. There seems to be 
no chance at all at the moment to get in- 
formation on the mantle Q from body tide 
gravity investigations . 

The delay angle of the potential bulge 
due to tidal deformation of the earth has 
been determined from the orbits of arti- 
ficial satellites to be 0.5° [Lambeck 
et al., 1974] . Lcimbeck et al. relate 

this delay angle to a mantle Q of 60 
which one obtains by using formula (3). 

The correct Q corresponding to this delay 
angle is 36, if the Gutenberg earth and 
constant Q values within the mantle are 
adopted. For the LMS equivalent model, 
we have calculated the global Q which 
corresponds to the delay angle of 0.5° to 
be 34. This shows that the delay angle of 
the tidal potential bulge is less sensi- 
tive to the geometry of the Q distribution 
within the earth than the phase delay of 
the tidal gravity at the deformed surface 
as measured by a gravimeter. The Q va- 
lues of 3 4 and 36 are much lov/er than the 
lowest limit of the upper mantle seismic 
Q. 

The usage of the exact theoretical 
relationship between the bulge of the 
tidal potential and the body's global Q 
may also be important for other planets 
as for instance for Mars. From observa- 
tions of the secular acceleration of 
the Mars satellite Phobos, Smith and 


simple formula (3), and found a global 
Q between 50 and 150. However, as one 
cannot neglect selfgravitation and the 
hydrostatic prestress for Mars, equation 
(3) is not applicable, and, therefore, 
the Q between 50 and 150 is probably too 
high for Mars, nrovided that the obser- 
ved phase angle of the tidal bulge is 
true . 

Similar considerations as above may 
also be important for the determination 
of the lunar global Q from observations 
of its physical librations [see Yoder, 
1978] . 

The Effect of Imoerfect Mantle Elasti- 
city on Loading Tides 

Unlike in the body tide case, also 
high degree harmonics of the load are 
important in the loading tide case. It 
is found that for the free oscillation 
Q model LMS, the loading tides of harmo- 
nic degrees slightly less than 100 are 
strongly effected by the low Q astheno- 
sphere, the global loading tide Qs being 
up to nearly 7 times smaller than the 
global body tide Q (see Fig. 2) . Corres- 
pondingly, the loading tide phase shifts 
due to imperfections in the elasticity 
of the mantle are by more than one order 
of magnitude higher than those of the 
body tide. Near ocean ridges and near 
subduction zones the M_ loading tide 
phase shifts may even Be as high as a 
several degrees for the disolacements as 
well as for gravity and tilt, provided the 
Maxwell constitutive law is valid. From 
the computation of phase shift Green's 
functions it is obvious that these maxi- 
mum phase shifts occur at about 80 to 
100 km distance from the load, this 
distance depending on the depth of the 
assumed low viscous asthenosnhere 
(see Fig. 3) . Our numerical results 
suggest that loading tide investigations 
could become an effective tool for 
studying the upper mantle viscosity in 
regions where viscosities lower than 
lol^ Poise may be expected. 
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Fig. 3. M 2 phase Green's functions for an 
ocean ridge model [ from Zschau, 1979a, b,]. 
They represent the phase shifts between the 
viscoelastic earth's responses to a varying 
point load, and the responses of the corres- 
ponding elastic earth. The high phase shifts 
at distances of about 100 km from the load 
are due to the low viscous asthenosphere. 

A !■ distance from the point load, PHI: 
phase shift. 
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